Binary Neutron Star Mergers: Dependence on the Nuclear Equation of State 
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We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear- 
theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore 
the dependence of the dynamical behavior of the binary neutron star merger and resulting gravita- 
tional waveforms on the EOS of the supernuclear-density matter. The numerical results show that 
the merger process and the first outcome are classified into three types; (i) a black hole is promptly 
formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is 
formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. 
For the EOS with which the maximum mass is larger than 2Mq , the lifetime of the HMNS is longer 
than 10 ms for a total mass mo = 2.7 Mq. A recent radio observation suggests that the maximum 
mass of spherical neutron stars is M max > 1.97 ± O.fMM© in one a level. This fact and our results 
support the possible existence of a HMNS soon after the onset of the merger for a typical binary 
neutron star with mo = 2.7M©. We also show that the torus mass surrounding the remnant black 
hole is correlated with the type of the merger process; the torus mass could be large, > O.IMq, in 
the case that a long-lived HMNS is formed. We also show that gravitational waves carry information 
of the merger process, the remnant, and the torus mass surrounding a black hole. 
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I. INTRODUCTION 



The coalescence of binary neutron stars is one of 
the most promising sources for kilometer-size laser- 
interferometric gravitational-wave detectors such as 
LIGO Q, VIRGO [2], and GEO600 0. A statisti- 
cal study based on the stellar evolution synthesis (e.g., 
Ref. Q) suggests that detection rate ~ 10 yr -1 may 
be achieved by advanced detectors such as advanced 
LIGO [|, advanced VIRGO @, and LCGT 7], which 
will be in operation in this decade. The merger of binary 
neutron stars will be a viable laboratory for studying 
supernuclear-density matter through gravitational-wave 
observations. For clarifying the nature as the sources of 
gravitational waves and for extracting their physical in- 
formation, theoretical templates of gravitational waves 
are needed. Because the merger is a highly dynamical 
process and proceeds in strong gravitational fields, nu- 
merical relativity is the unique way for deriving the tem- 
plates of gravitational waves. 

The merger of binary neutron stars also has been pro- 
posed as a likely candidate for the central engine of 
short 7-ray bursts (GRBs) [1, @. The observations have 
shown that the central engine supplies a large amount 
of energy > 10 48 ergs in a short time scale < 2 s [lol ]. 
According to a standard scenario based on the merger 
hypothesis, a stellar-mass black hole surrounded by a 
hot massive torus should be formed after the merger. 
Possible relevant processes to extract the energy of this 
black hole-torus system for launching a relativistic jet are 
neutrino- antineutrino pair annihilation and/or magneti- 
cally driven mechanisms. Recent numerical studies (e.g., 
Ref. [ll|) suggest that if the torus has a mass > O.IMq, 
it could supply the required energy of short GRBs by 
the neutrino mechanism. The amount of the mass of the 



remnant torus depends on the dynamical behavior of the 
merger process of binary neutron stars. Thus, the issue 
is to clarify how the formation process of a massive torus 
depends on the dynamics of the merger, on the equation 
of state (EOS) of neutron stars, and on parameters of 
the binary such as total mass and mass ratio, for under- 
standing the formation mechanism of the central engine 
of short GRBs. 

The supernuclear-density EOS plays a key role for de- 
termining the merger process of binary neutron stars. 
For example, for a given mass, the outcome of the 
merger depends strongly on the EOS: For soft EOSs, 
the merger results in prompt formation of a black hole. 
On the other hand, for stiff EOSs, a hypermassive neu- 
tron star (HMNS) is formed [Ea - tl4l |. However, the actual 
supernuclear-density EOS is still unknown because of our 
poor knowledge about the properties of the matter above 
the nuclear density. 

In this paper, we extend the previous works for a long- 
term simulation of binary neutron stars (e.g., Refs. [15l - 
[j~7l |) from the following motivation. Recently, a piecewise- 
polytropic EOS for the cold EOS, based on the nuclear 
theoretical calculations, was pro posed by Read et ai, 
and by Ozel and Psaltis [TH Il9j|. With only four pa- 
rameters, the piecewise-polytropic EOS can approximate 
any candidate of the EOS of supernuclear-density mat- 
ter. By using this EOS, we can systematically study the 
effects of the possible EOSs on phenomena related to the 
supernuclear-density matter, such as the merger of bi- 
nary neutron stars. 

In this work, we report our latest numerical-simulation 
results for the merger of equal-mass binary neutron stars, 
for which the total masses are in the range of 2.7 Mq- 
3.0M Q . To systematically study the effects of the EOS 
on the merger, we use 6 different stiff EOSs which are de- 
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scribed in Ref. 18]. In the present simulation, we follow 
the inspiral motion for 5-7 orbits and the merger process 
up to formation of a stationary black hole or HMNS. Af- 
ter a black hole is formed, we study the dependence of 
the properties of the resulting torus, such as its mass, on 
the EOS and the total mass of the binaries. We also ex- 
tract gravitational waves and show the dependence of the 
gravitational waveforms and their spectra on the EOS. 

The paper is organized as follows. In Sec. II, we sum- 
marize the formulation and numerical schemes employed 
in our numerical code SACRA, and the EOS employed in 
this study. In Sec. III. A, we describe our simulation re- 
sults for the typical cases of the merger. We define three 
types of the merger process, which clearly classify the 
dependence of the dynamical behavior of the merger on 
the EOSs. In Sec. III.B, we summarize the character- 
istic features of gravitational waveforms for each type. 
Section IV is devoted to a summary. Throughout this 
paper we use the geometrical units of c = G = 1 where c 
and G are the speed of light and gravitational constant, 
respectively. 

II. FORMULATION 

A. Numerical method 

We follow the late inspiral and merger phases of binary 
neutron stars using a numerical-relativity code, called 
SACRA, described in Ref. [2(j. SACRA employs a mov- 
ing puncture version of the Baumgarte-Shapiro-Shibata- 
Nakamura formalism 21 23] to solve Einstein's evolu- 
tion equation without imposing any symmetry. In SACRA, 
we evolve a conformal factor W = 7" 1 / 6 , the conformal 
three metric 7^ = 7~ 1 / 3 7ij, the trace of extrinsic cur- 
vature K, the conformal trace-free extrinsic curvature 
Aij = (Kij — K^/ij/3), and an auxiliary variable 

V = — 9j7 y . Here 7^ is the three metric, Kij is the 
extrinsic curvature, and 7 = det(7ij). In the numerical 
simulation, a fourth-order finite differencing scheme in 
space and time is used implementing an adaptive mesh 
refinement (AMR) algorithm (at refinement boundaries, 
a second-order interpolation scheme is partly adopted). 
The advection terms such as ^di^jk are evaluated by a 
fourth-order non-centered finite difference 24] . A fourth- 
order Runge-Kutta method is employed for the time evo- 
lution. 

Following Ref. [24| , we adopt a moving-puncture gauge 
condition as 

(d t - pdj)F = 0.75^', (1) 
(d t - p j d j )B i = (d t - &dj)f* - i h B\ (2) 

where B l is an auxiliary variable and rj s is an arbitrary 
constant. In the present paper, we set 77,5 ~ 3/M. Here, 
M denotes the mass for each neutron star in isolation. 

For the hydrodynamics, we employ a high-resolution 
central scheme by Kurganov and Tadmor [251 ] with a 



third-order piecewise parabolic interpolation and with a 
steep min-mod limiter. 

In SACRA, an AMR algorithm is adopted (see Ref. [20| 
for details). In the present work, we prepare seven re- 
finement levels both to accurately resolve the structure 
of two neutron stars and to extract gravitational waves 
in a local wave zone. In our simulations, two sets of 
four finer domains comoving with the neutron stars cover 
the region in their vicinity. The other three coarser 
domains cover both neutron stars by a wider domain 
with their origins being fixed approximately at the cen- 
ter of mass of the binary. Each refinement domain con- 
sists of the uniform, vertex-centered Cartesian grids with 
(2N + 1, 2N + 1, N + 1) grid points for (x, y, z) with the 
equatorial plane symmetry at z — imposed. The half of 
the edge length of the largest domain (i.e., the distance 
from the origin to outer boundaries along each axis) is de- 
noted by L which is chosen to be > Ao, where Ao = 7r/f2o 
is the initial wavelength of gravitational waves. The grid 
spacing for each domain is then hi = L/(2 l N), where 
1 = 06. In this work, we typically choose TV = 60. 
With this grid resolution, the semi-major diameter of 
each neutron star is covered by about 100 grid points 
(cf. Table HI)) . In addition, we performed numerical sim- 
ulations with lower grid resolutions, N = 36, 42, and 50, 
to check the convergence of the numerical results. The 
property of the convergence is essentially the same as in 

Refs. dam. 

B. Models of Equation of State 

The parameterized piecewise-polytropic EOS [HI, [l9| 
is useful to systematically study the dependence of the 
dynamical behavior of the merger on the EOS of the 
supernuclear-density matter. In this work, we employ 
a parameterized piecewise-polytropic EOS proposed by 
Read et al. [HI. This EOS is written in terms of four 
segments of polytropes 

P =K iP Ti (3) 
( for Pi < p < pi+i, < i < 3), 

where p is the rest-mass density, P is the pressure, Ki 
is the polytropic constant, and Ti is the adiabatic index. 
We refer to the pressure in the form of Eq. (3) as the 
cold-part pressure, P C oid- At each boundary of the piece- 
wise polytropes, p = pi, the pressure is required to be 
continuous, i.e., K.ip\' = Ki + ip^' +1 . Read et al. de- 
termine these parameters in the following manner [l8j . 
First, they fix the EOS of the crust as Tq = 1.357 and 
Kq = 3.594 x 10 13 in the cgs unit. Then they determine 
p 2 = 1.85p nuc i and p 3 = 3.70p nuc i where p nuc \ = 2.7x 10 14 
g/cm 3 is the nuclear saturation density. With this prepa- 
ration, they choose the following four parameters as a 
set of free parameters: {Pi, Ti, T 2 , T 3 }. Here Pi is the 
pressure at p — p%, and from this, K\ and Ki are deter- 
mined by Ki = Pi/p\ x and ifj+i = Kip\ l r,+1 . There- 
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TABLE I: Parameters of the piecewise-polytropic EOS, the maximum mass of spherical neutron stars, M max , and the radius 
of a spherical neutron star of mass M = 1.4M Q , R1.4, for each EOS. Composition means strongly interacting components 
(n=neutron, p=proton, H=hyperon, Q=quark, 7r°=pion) and APR4, SLy, H3, H4, and ALF2 include leptonic components. 



EOS 


log Pi ( dy ne/ cm 2 ) 


Ti 


r 2 


r 3 


A/max (M ) 


i?i.4(km) 


Approach 


composition 


APR4 


34.269 


2.830 


3.445 


3.348 


2.213 


11.428 


Variational-method 


np 


SLy 


34.348 


3.005 


2.988 


2.851 


2.049 


11.736 


Effective-one-body potential 


np 


H3 


34.646 


2.787 


1.951 


1.901 


1.788 


13.840 


Relativistic mean field 


npH 


H4 


34.669 


2.909 


2.246 2.144 


2.032 


13.759 


Relativistic mean field 


npH 


ALF2 


34.055 


4.070 


2.411 


1.890 


2.086 


13.188 


APR+ Quark matter 


npQ 


PS 


34.671 


2.216 


1.640 


2.365 


1.755 


15.472 


Pion condensation 


nvr 



fore the EOS is specified by choosing the four parameters 

{Px.rx.ra.ra}. 

In this paper, we adopt 6 models of piecewise- 
polytropic EOS which describe the following EOSs based 
on nuclear theoretical calculations. 

1. APR4: derived by a variational-method with the 
AV18 2-body potential, the UIX 3-body potential, 
and relativistic boost corrections (see Ref. [26j): 

2. SLy: derived by using an effective potential ap- 
proach of the Skyrme type (see Ref. [271]); 

3. H3: derived by a relativistic mean-field approach 
including hyperons. The incompressibility, the ef- 
fective mass, and the nucleon-meson coupling are 
chosen to be K = 300 MeV, m*/m„ = 0.7, and 
x a = 0.6. Here m n is the nucleon-mass (see 
Refs. dHH); 

4. H4: the same as H3 but for x a — 0.72 (see Refs. [28], 
HI); 

5. ALF2: a hybrid EOS which describes nuclear mat- 
ter for a low density and color-flavor-locked quark 
matter for a high density. The transition density 
and the interaction parameter are chosen to be 
Pc = 3p nuc i and c = 0.3 (see Ref. [3(ij): 

6. PS: derived by using a potential approach. This 
EOS describes a neutron matter with pion conden- 
sation (see Ref. |3lj|). 

TableUlists the parameters of piecewise-polytropic EOSs 
employed in this work. These EOSs are relatively stiff, 
and hence, the maximum mass of spherical neutron stars 
is larger than 1.75M . The choice of these EOSs is mo- 
tivated by the recent discovery of a heavy neutron star 
with mass 1.97 ± 0.04A/ Q (one a error) [32J. This value 
has become the new standard for the minimum value of 
the neutron star maximum mass }33l l34| . 

The thermal pressure should be taken into account 
for numerical simulations, because matter in the merged 
neutron stars becomes hot with temperature to T > 10 
MeV due to the shock heating at the merger (e.g., 
Refs. [H, In this case, the thermal energy is not 



negligible. To approximately include the thermal pres- 
sure, we employ the EOS which is described by 

P(p,s)=P cold (p) + P th (p,e), (4) 

where e is the specific internal energy, P co id is the pres- 
sure determined by the piecewise-polytropic EOS, and 
Pth is the thermal part of the pressure which is given by 

Ph = (T th - 1) (e - £coid)p. (5) 

Here e co id is determined from P co id by the first law of 
thermodynamics [l8| . In our simulations, we focus only 
on the case that the shock heating efficiency is relatively 
low, i.e., r t h = 1.357. To study the effect of thermal 
pressure on the outcome of the merger, we also employed 
r t h = 1.5, 1.7, and 1.8 for a few simulations with lower 
grid resolutions. We found that the outcome formed soon 
after the merger (a black hole or HMNS) depends very 
weakly on the value of r t h (see also Ref. [13j ) , although 
the long-term evolution process of a HMNS depends on 

it m. 

C. Initial data 

We prepare binary neutron stars in quasiequilibrium 
states for the initial condition of numerical simulations. 
To track more than 5 quasicircular orbits with small ec- 
centricity for deriving accurate gravitational waveforms 
in the late inspiral and merger phases, orbital separation 
of the initial configuration is chosen to be large enough 
that the time scale of gravitational radiation reaction is 
much longer than the orbital period of the binary. 

The formulation and methods for a solution of Ein- 
stein's constraint equation and equations of hydrostat- 
ics are the same as those adopted in our previous 
works HE US EH, except for the choice of EOSs. We 
assume the conformal flatness of the three metric 7^ = 
ip fij, the maximal slicing condition K — 0, and their 
preservation in time. Here, tp is a conformal factor and 
fij is the fiat spatial metric. The piecewise-polytropic 
EOS described in the previous subsection is adopted to 
model the neutron star matter because the fluid inside 
the neutron stars in the late inspiral phase are believed to 
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TABLE II: Key parameters for the initial models adopted in the numerical simulation, mo is the sum of the ADM masses of 
two neutron stars in isolation (2M); M ADM and Jq" um are the ADM mass and angular momentum of the system, respectively; 
M* is the baryon rest mass; Qo is the angular velocity. We also show the setup of the grid structure for the computation with 
our AMR algorithm. Ax = he = L/(2 6 N) is the grid spacing at the finest resolution domain with L being the location of the 
outer boundaries for each axis. -Rdiam denotes the number of the grid points assigned inside the semimajor diameter of the 
neutron stars. Ao is the wevelength of gravitational waves of the initial configuration. In the last two columns, we show the 
simulation results for the rest mass of the torus M t0 rus and the type of the merger process for each model. The rest mass of 
the torus surrounding the black hole is determined at 1 ms after the black hole formation. For APR4-27, a HMNS with the 
lifetime ^> 10 ms is formed. For H4-30, the type of the merger process is ambiguous because a black hole is formed at only 1.5 
ms after the onset of the merger. 



Model 


TTIq 


M„ ADM 


tAUM 

"0 


M* 


Win Qn 


At J TTIq Rd 


iam / * — 1 


L/Ao 


Aftnr,, a /Mn 


tvDe 


APR4-27 


2.7 


2.67 


7.16 


3.00 


0.0221 


0.043 


99 


1.16 




in 


APR4-28 


2.8 


2.77 


7.70 


3.12 


0.0221 


0.041 


102 


1.11 


0.003 


i 


APR4-29 


2.9 


2.87 


8.26 


3.26 


0.0221 


0.039 


102 


1.05 


<0.001 


i 


SLy-27 


2.7 


2.67 


7.16 


2.98 


0.0221 


0.045 


101 


1.21 


0.02 


ii 


SLy-28 


2.8 


2.77 


7.70 


3.12 


0.0221 


0.043 


102 


1.15 


< 0.001 


i 


H3-27 


2.7 


2.68 


7.39 


2.94 


0.0221 


0.056 


102 


1.50 


0.05 


ii 


H3-29 


2.9 


2.87 


8.27 


3.18 


0.0221 


0.050 


103 


1.34 


0.01 


i 


H4-27 


2.7 


2.68 


7.39 


2.94 


0.0221 


0.056 


103 


1.50 


0.18 


in 


H4-29 


2.9 


2.87 


8.27 


3.18 


0.0221 


0.051 


101 


1.37 


0.02 


ii 


H4-30 


3.0 


2.97 


8.85 


3.30 


0.025 


0.048 


102 


1.49 


0.01 


I or 11 


ALF2-27 


2.7 


2.67 


7.17 


2.98 


0.0221 


0.049 


102 


1.32 


0.16 


III 


ALF2-29 


2.9 


2.87 


8.51 


3.22 


0.0221 


0.045 


102 


1.22 


0.02 


II 


ALF2-30 


3.0 


2.97 


8.85 


3.34 


0.0221 


0.043 


102 


1.32 


0.003 


I 


PS-27 


2.7 


2.68 


7.57 


2.88 


0.020 


0.073 


92 


1.60 


0.04 


III 


PS-29 


2.9 


2.88 


8.73 


3.12 


0.020 


0.065 


92 


1.48 


0.02 


II 


PS-30 


3.0 


2.97 


8.85 


3.24 


0.025 


0.056 


102 


1.71 


0.01 


I 



be well approximated by a cold, zero-temperature mat- 
ter. The neutron stars are assumed to have an irrota- 
tional velocity field, which is believed to be an astrophys- 
ically realistic configuration [3^, [40] . Numerical compu- 
tations are performed using the spectral-method library, 
LORENE [4l| . The details of the numerical methods and 
the analysis of quasiequilibrium states are summarized in 
Ref . (see also Ref . [38[ ) . We note that the virial error 
of the quasiequilibrium, which we define as the relative 
difference between the ADM and Komar masses, is al- 
ways smaller than 10 -4 for our numerical solutions. 

D. Gravitational wave extraction and waveforms 

Gravitational waves are extracted by calculating the 
complex Weyl scalar $4. using the same procedure as in 
Ref. [2(| ■ Gravitational waveforms are calculated by 

h+(t)-ih x (t) = - lim / dt' [ dt"V A (t" ». (6) 

r^ooj J 

Here we omit arguments 9 and cf>. We evaluate ^4 at 
a finite coordinate radius r = 400M© ~ 590 km. In 
the standard spherical coordinate (r, 8, </>), ^4 can be ex- 
panded in the form 

f 4 (i, r, 6,<f>) = J2 *4 m (^ r)_ 2 Y lm (9, 0), (7) 

in 1 



where -2Yi m are spin- weighted spherical harmonics of 
weight —2 and W™ are expansion coefficients defined 
by this equation. In this work, we focus only on the 
(I, |m|) = (2,2) mode. 

We evaluate the amplitude of the Fourier spectrum of 
gravitational waves, 

Hn = / M/)l 't few ' 2 . w 

where / is the frequency, and h + (f) and h x (f) are the 
Fourier transformation of the plus and cross modes of 
gravitational waves observed along the z axis. The ef- 
fective amplitude of gravitational waves for a given fre- 
quency is defined by, 

m/) = fkf)- (9) 

Note that this is the most optimistic value for the effec- 
tive amplitude. Actually, the amplitude of gravitational 
waves depends on an angle locating the source in the sky 
and on an angle specifying the orientation of orbital plane 
of the binary neutron star. The angular average of the 
effective amplitude is approximately ~ 0.4ft, c ff. 

E. Mass, linear momenta and angular momenta 

We monitor the ADM mass A^adm , the linear momen- 
tum Pi, and the angular momenta Ji during the evolu- 
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tion. These parameters are defined by the integrals on 
two surfaces of a coordinate radius r = 400, 300, and 
240M Q , 

M ADM (r) = ^ J Vl7 ij l kl (lik,i ~ lij,k)dS u (10) 
Pi{r) = ^fSr&t-KTbdSj, (11) 
Mr) = ^e m J^x\K^ k -K^ k )dS„ (12) 

where dSi is the surface element and e^fc is the Levi- 
Civita symbol. Then, we extrapolate these quantities for 
r — > oo to obtain the asymptotic value. 

We also monitor the total baryon rest mass 



because a typical total mass of binary neutron stars is 



2.7 M P 



© 



-gd 3 x, 



(13) 



where u* is the time-component of the four velocity, and 
g is the determinant of the space-time metric. After the 
black hole formation, we calculate the torus mass defined 
by 

Mtorus = / pu t y^gd 3 X, (14) 

Jr>r A n 

where tah is the coordinate radius of the apparent hori- 
zon. 



III. NUMERICAL RESULTS 

Table HTl lists the numerical models adopted in this pa- 
per. The simulations were performed from 5-7 orbits 
before the onset of the merger to 3 ms after the formation 
of a black hole or to the phase in which a HMNS relaxes 
to a quasi-stationary state. 



A. Dynamics 

The merger dynamics are determined primarily by 
three forces: gravity, pressure, and centrifugal forces. If 
its gravity is stronger than the other forces, the merged 
neutron star collapses to a black hole soon after the on- 
set of the merger. On the other hand, if the other forces 
could overcome the gravity in a compact state, the core 
of the merged neutron star bounces and a rapidly ro- 
tating and oscillating HMNS is formed. This dynamical 
behavior depends on the EOS and on the total mass of 
the binary neutron star, tuq. In the following, we show 
the results of our numerical simulations focusing on the 
dynamical behavior of the merger. 



1. Classification of the merger process 

It is natural to expect that a black hole is eventu- 
ally formed after a merger of any binary neutron star, 



42j . which is likely to exceed the maximum 
mass of spherical neutron stars. Although binary neu- 
tron stars have such a large total mass, a HMNS is often 
formed as a transient outcome of the merger, which is 
supported by the strong centrifugal force caused by rapid 
and differential rotation as well as by thermal pressure. 
Thus, one can classify the merger process and the result- 
ing remnant of binary neutron stars into three types (see 
Table II): 

• type I: A black hole is promptly formed; 

• type II: A short-lived HMNS is formed 
(t h < 5 ms); 

• type III: A long-lived HMNS is formed 
(th > 5 ms). 

Here th is a lifetime of a HMNS. We note that a massive 
and stable rigidly rotating neutron star may be formed 
if the total mass is only slightly larger than the maxi- 
mum mass of spherical neutron stars or the contribution 
of the thermal pressure is significant. However, we do 
not consider this fourth possibility in this paper (but see 
Ref. [li|). 

Figure [T] shows the snapshots of the density colormap 
for three types of the merger process. Since the dynamics 
of the binary neutron stars in the inspiral phase is similar 
among three types, we focus only on the dynamics after 
the onset of the merger in the following. 

Top panels; APR4-29 (type I). Soon after the onset of 
the merger (Fig.[TJ top left,), the merged object collapses 
promptly to a black hole and tiny materials remain out- 
side the black hole (Fig. (TJ top center). The black filled 
circle denotes the inner region of the apparent horizon. 
Note that the spiral arms are formed and the materials 
in their outer region obtain angular momentum from the 
materials in the inner region by gravitational torques re- 
sulting from the non-axisymmetric structure. At ~ 5 ms 
after the onset of the merger, a quasi-stationary torus is 
formed with the maximum density p max ~ 10 11 ' 5 g/cm 3 , 
and spreads to about 10 km from its center (Fig. [TJ top 
right). However the torus mass is small < 1O _3 M . 

Middle panels; H3-27 (type II). After the onset of 
the merger, the merged core bounces due to strong cen- 
trifugal and pressure forces. The resulting HMNS has a 
double-core structure in which they rotate around each 
other (Fig.[TJ middle center). At ~ 5 ms after the onset 
of the merger, the HMNS collapses to a black hole (Fig.Q] 
middle right), because its angular momentum decreases 
due to the emission of gravitational waves. We find that 
the collapse of the HMNS occurs during the phase in 
which the HMNS has a non-axisymmetric shape. The re- 
sulting torus around the black hole (Fig. Q] middle right) 
spreads to about 20 km. This torus has the maximum 
density p ma x ~ 10 12 5 g/cm 3 . In this case, the torus mass 
is ~ 0.05-M©. The resulting black hole-torus system is a 
candidate for the central engine of short GRBs. 
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FIG. 1: Colormap of the density, logp (g/cm 3 ). Top, middle, and bottom rows show the snapshots for APR4-29, H3-27, and 
H4-27, respectively. The black filled circle denotes the region inside apparent horizon. Note that the density range of the color 
bar for APR4-29 is different from the other models. 



For H4-29, ALF2-29, and PS-29, a HMNS of double- 
core structure is also formed while for SLy-27 a HMNS of 
ellipsoidal shape is formed. Irrespective of the configura- 
tion of the HMNS, however, the mass of remnant torus 
formed after the black hole formation is ~ 0.02Af Q and 
much less than that for H3-27. This suggests that the 
torus mass may be suppressed for a massive system, for 
which the total mass is close to M cr it; see Eq. (15) for 
the definition of M cr i t . 



Bottom panels; H4-27 (type III). After the onset of 
the merger, a HMNS, which has a double-core structure, 
is formed as in the case of H3-27 (Fig. [TJ bottom left). 
Because the lifetime of the HMNS is sufficiently long, 
a large amount of angular momentum is dissipated by 
gravitational waves. As a result, the HMNS approaches 
an axisymmetric ellipsoidal shape (Fig. [T] bottom cen- 
ter). After the significant loss of the angular momen- 
tum, it collapses eventually to a black hole (Fig. [1] bot- 
tom right). The resulting torus surrounding the formed 
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FIG. 2: The evolution of the maximum baryon rest-mass den- 
sity, Pmax, for three models. The solid, dashed, dash-dotted, 
and dotted curves denote the results for models APR4-29 
(type I), H3-27 (type II), H4-27 (type III), and APR4-27 (type 
III), respectively. 



black hole has the maximum density p max ~ 10 13 g/cm 3 
and spreads widely to r ~ 25 km. In this case, the torus 
mass is ~ 0.18Af Q . The resulting black hole-torus system 
may be regarded as a promising candidate for the central 
engine of short GRBs. 

APR4-27, ALF2-27, and PS-27 also show the merger 
process of type III. For ALF2-27, the mass of the rem- 
nant torus is larger than 0.1M Q as for H4-27. This sug- 
gests that with the stiff EOSs for which the maximum 
mass of spherical neutron stars is larger than 2M Q , a 
massive torus could be the outcome when the total mass 
of the binary is « 2.7M . For PS-27 with which the 
maximum mass of spherical neutron stars is 1.75M©, de- 
spite type III, the mass of the remnant torus is much less 
(~ 0.04M©). This indicates that for the EOS with small 
maximum mass, the torus mass may be suppressed. More 
detailed reasons that could cause the difference in the 
remnant torus mass will be discussed in Sec. III. A. 3. For 
APR4-27, a quasi-stationary HMNS is formed. Because 
the degree of differential rotation is still high and the 
mass is much higher than the maximum mass of spherical 
neutron stars, it will eventually collapse in the presence 
of the magnetic field or viscosity in reality (see Ref . T3!] ) . 

Figure [2] plots the evolution of the maximum baryon 
rest-mass density, p max , for models (types) APR4-29 
(type I), H3-27 (type II), H4-27 (type III), and APR4- 
27 (type III). For APR4-29, for which a black hole is 
formed promptly, p max increases monotonically after the 
onset of the merger. For H3-27 and H4-27 for which a 
HMNS is formed transiently, p max oscillates and then in- 
creases until the collapse of the HMNS sets in. After the 
collapse of the HMNS, p max increases monotonically. For 
APR4-27, p max also oscillates and increases soon after the 
formation of the HMNS. However it eventually reaches a 



relaxed value, implying that a nearly stationary HMNS 
is the outcome. 

Figure [3] summarizes the type of the merger process for 
each EOS and for each mass. We find that a HMNS is 
likely to be formed for the merger with the EOS which 
provides a small compactness, M/R, for each neutron 
star. The reason is as follows. In the case that the com- 
pactness of each neutron star is small, two neutron stars 
merge at a relatively large orbital separation. As a re- 
sult, the merged neutron stars have large angular momen- 
tum at the onset of the merger, which helps escaping the 
prompt collapse to a black hole. The long-lived HMNS is 
also likely to be formed for a total mass tuq > 2.7M Q with 
the EOS which has the maximum mass exceeds 2M@, 
such as APR4, H4, and ALF2. 

Note that it is practically impossible to precisely de- 
termine the lifetime of the HMNS by the numerical sim- 
ulation, because the HMNS just before the collapse is 
marginally stable and its dynamics depends strongly on 
a small perturbation and thus on the grid resolution. The 
lifetime also depends on the treatment of the thermal ef- 
fects (3a |. which are determined by the value of T t h in 
our simulations. However we find that our classification 
of the merger process depends very weakly on the grid 
resolution and the value of F t h- 



2. The Critical mass 

We define a critical mass, M cr j t , of binary neutron stars 
for each EOS as follows: If the total mass of a binary 
neutron star, mo, exceeds M cr it, a black hole is formed 
promptly after the onset of the merger. 

We write the value of M cr it in terms of M max as 



M cr it = kM n 



(15) 



where A; is a constant which depends on the EOS. Ta- 
ble |HT] shows the maximum mass, the critical mass, and 
their ratio, k, for each EOS together. Here we find that 
k is in the range, 



1.3 < k < 1.7, 



(16) 



where k ~ 1.3 for APR4 (R 1A = 11.428 km) and SLy 
(Ria = 11.736 km) and k ~ 1.7 for PS (R 1A = 15.472 
km). We find the correlation between k and Ria, in 
which k is approximately an increasing function of Ria- 
Note that the results for APR4 and SLy agree with those 
in the previous papers [H, [lj[ • 

Recent observation suggests M max > 1.97±O.O4M in 
one a level [Hj]. This fact and our results suggest that 
a HMNS is likely to be formed transiently soon after the 
onset of the merger for a binary neutron star merger with 
the typical total mass, toq — 2.7M Q [42| . 
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FIG. 3: Type of the merger process and the remnants for each model. The vertical axis shows the total mass of two neutron 
stars. The horizontal axis shows the EOSs together with their radii for M — 1.4M©, R\a km. 



TABLE III: The maximum mass, M max , the critical mass, M cr it, and their ratio k for each EOS. 





APR4 


SLy 


ALF2 


H3 H4 


PS 


k 


1.3 


1.3 


1.4 


1.6 1.5 


1.7 


M max /M Q 


2.213 


2.049 


2.086 


1.788 2.032 


1.755 


M crit /M Q 


~ 2.8 


~ 2.8 


~ 3.0 


~ 2.9 > 3.0 


~ 3.0 



3. Final states: Black hole and Torus 

As summarized in Sec. III.A.l, after the merger of a 
binary neutron star, a black hole surrounded by a torus 
is eventually formed (see Fig. [T]). Our interest here is to 
study how the amount of the torus mass depends on the 
type of the merger process. Figure |4] shows the evolution 
of the torus mass for models (types) APR4-29 (type I) , 
H3-27 (type II), and H4-27 (type III). Here the time at 
the black hole formation is set to be t = 0. Table ILT1 shows 
the results for the torus mass as well as the type of the 
merger process for each model. We find that the torus 
mass is correlated to the type of the merger process as 
follows, 

M t or us 

< 0.01M Q for type I, (17) 
0.02A/ Q < Aftorus < 0.05M© for type II, (18) 
0.04M Q < M torus < 0.18M© for type III. (19) 



Here the torus mass is evaluated at about 1 ms after 
the black hole formation. Thus, the torus mass is larger 
in the case that a HMNS is formed than that a black 
hole is formed promptly. When a long-lived HMNS is 
formed, in particular, the resulting torus mass could be 
Mtorus > 0.1M© for many EOSs. This feature can be un- 
derstood as follows. In the HMNS, its envelope spins up 
because angular momentum is transported from inner to 
outer parts of the HMNS by gravitational torques caused 
by the non-axisymmetric structure. As a result, a part 
of the matter in the envelope of the HMNS, which has 
sufficient angular momentum, does not fall into the black 
hole at the collapse of the HMNS and it constitutes the 
torus. However, we note that the efficiency of the angu- 
lar momentum transport depends on the density profile 
and the degree of non-axisymmetry of the HMNS, which 
depend on the EOS and the total mass. For PS-27, in- 
deed, the resulting torus mass is only s» 0.04A/© which 
is much less than those for APR4-27, H4-27, and ALF2- 
27. Thus, for a special EOS which leads to a large radius 
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and a small maximum mass such as R\a — 15.5 km and 
M-caax — 1-76M Q , the torus mass may be small even if 
the merger process is type III. 




FIG. 4: The evolution of the torus mass, M tor us, for three 
models. The solid, dashed, and dash-dotted curves denote 
the results for models APR4-29 (type I), H3-27 (type II), and 
H4-27 (type III), respectively. The time at the black hole 
formation is set to be t = 0. 

Note that, in this study, we performed numerical sim- 
ulations only for the equal-mass systems. For unequal- 
mass systems, a massive torus may be formed even for the 
case that the merged neutron stars collapse promptly to 
a black hole. This is because the heavier star can disrupt 
the less massive companion by tidal forces, in particu- 
lar, for high mass ratio and subsequent angular momen- 
tum transport in the spiral arms formed from the tidally 
disrupted neutron star enhances the torus formation, as 
indicated in Refs. H-El. 



B. Gravitational Waves 

Gravitational waves are emitted during the merger un- 
til a stationary black hole is formed. The gravitational 
waveforms reflect the dynamical behavior of the merger 
process. In the following, we classify the gravitational 
waveforms and their spectra into three types in the same 
way as the merger process, and discuss their features. 



1. Gravitational waveforms 

The gravitational waveforms for models APR4-27, H3- 
27, H4-27, and PS-27 together with post-Newtonian 
waveforms, calculated by the so-called Taylor T4 for- 
mula 47j, are shown in Figs. [5] (la) - (4a). For more 
massive models APR4-29 and H4-29 the gravitational 
waveforms are shown in Figs. [5] (5a) and (6a). The grav- 
itational waveform for the inspiral phase is similar for all 



the models. In the late inspiral stage, the gravitational- 
wave phase starts deviating from the post- Newtonian one 
at 3 - 5 ms before the onset of the merger. This is because 
the effects of the finite size of two neutron stars play an 
important role for their orbital motion. A detailed anal- 
ysis for the finite-size effect will be published in a future 
paper. In the following, we focus on the gravitational 
waveform of the merger and the ringdown phases. 

Type I. The gravitational waveform for APR4-29 is 
shown in Fig. [5] (5a). The amplitude of gravitational 
waves increases gradually in the inspiral phase until the 
merger sets in at t ~ 22 ms. Soon after the onset of the 
merger, ringdown gravitational waves are emitted by the 
oscillating black hole for ~ 1 ms. Then the amplitude 
approaches zero because the black hole becomes station- 
ary. 

Type II. A short-lived HMNS is formed after the onset 
of the merger. Then the oscillating and rotating HMNS 
emits quasiperiodic gravitational waves. The gravita- 
tional waveform of H3-27 is shown in Fig. [5] (2a). In this 
case, two neutron stars merge at t ~ 18 ms, at which 
the amplitude of gravitational waves is small transiently. 
This implies that the merged object has a nearly ax- 
isymmetric ellipsoidal shape at the moment. After this, 
the core bounces and the HMNS is formed, which has a 
double-core structure. Then quasiperiodic gravitational 
waves with a high amplitude are emitted. At t ~ 23 
ms, the amplitude of gravitational waves damps sud- 
denly. This is because the HMNS collapses to a black 
hole before the HMNS becomes an axisymmetric ellip- 
soidal shape. A similar waveform is also seen for H4-29 
in Fig. M (6b). 

Type III. The gravitational waveform of H4-27 is shown 
in Fig. [5] (3a). In this case, two neutron stars merge 
at t ~ 18 ms, at which the amplitude of gravitational 
waves is small transiently as in H3-27. After this, a 
non-axisymmetric HMNS of the double-core structure is 
formed and it emits quasiperiodic gravitational waves. 
For t > 30 ms, gravitational waves of small amplitude is 
emitted quasi-stationarily until t ~ 37 ms. During this 
phase, the HMNS has a nearly axisymmetric ellipsoidal 
shape. At t ~ 37 ms, the HMNS collapses to a black hole 
and the amplitude of gravitational waves damps eventu- 
ally. 

For APR4-27 (see Fig. (la)), an ellipsoidal HMNS is 
formed after the onset of the merger, and thus, quasiperi- 
odic gravitational waves are emitted for the first ~ 10 
ms after the formation of the HMNS. The HMNS loses 
angular momentum due to the gravitational radiation re- 
action and the ellipticity gradually decreases, resulting in 
the decrease of the gravitational-wave amplitude. Thus, 
the gravitational waveform is similar to that for H4-27. 
However, the lifetime of the HMNS for APR4-27 is much 
longer than 15 ms. Thus we did not follow the collapse 
of the HMNS for APR4-27. 

For PS-27 (see Fig. [5] (4a)), the gravitational wave- 
form is different from those of H4-27 and APR4-27. In 
this case, it is similar to those for H3-27 and H4-29: 
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FIG. 5: Gravitational waveforms and their spectra. The solid and dashed curves in the left panels denote the waveforms 
calculated by the simulation and Taylor T4 formula, respectively. The solid and dashed curves in the right panels denote the 
spectra calculated by the simulation, and spectrum calculated by Taylor T4 formula, respectively, at a hypothetical source 
distance of 100 Mpc. The effective amplitude for the most optimistic direction of the source is shown. Here the noise levels of 
advanced LIGO (Optimal NSNS version), LCGT (Broadband version), and Einstein Telescope are shown together, (la) and 
(lb) for APR4-27 (type III), (2a) and (2b) for H3-27 (type II), (3a) and (3b) for H4-27 (type III), (4a) and (4b) for PS-27 
(type III). 
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FIG. 6: The same as Fig. |H but for m = 2.9M . (5a) and (5b) for APR4-29 (type I), (6a) and (6b) for H4-29 (type II). 



A HMNS of double-core structure is formed after the 
onset of the merger, and emits quasiperiodic gravita- 
tional waves. After a substantial emission of gravita- 
tional waves, the HMNS collapses to a black hole before 
it becomes a nearly axisymmetric ellipsoid. 



2. Spectra 

The right panels of Figs. [5] and [6] show the spectra of 
gravitational waves for models APR4-27, H3-27, H4-27, 
PS-27, APR4-29, and H4-29 together with the spectrum 
calculated by the post-Newtonian approximation (Tay- 
lor T4), and the sensitivity curves of advanced LIGO 
(Optimal NSNS version) gjj, of LCGTjBroadband ver- 
sion) [7j, and of Einstein telescope |49j |. Here we as- 
sume that gravitational waves are observed at a distance 
of 100 Mpc from the sources along the most optimistic 
source direction. The common feature for each type is 
that the effective amplitude, h e g, decreases with the in- 
crease of / until / — 1 kHz, in the inspiral phase. The 
spectrum shape above / ~ 7 kHz is also qualitatively 
the same irrespective of the model, which is caused by 
the quasi-normal-mode oscillations of the remnant black 
hole. However, the spectrum shape between 1 kHz and 
7 kHz, for which gravitational waves are emitted in the 
merger phase, depends strongly on the EOS and on the 
total mass as summarized in the following: 



Type I. Because no HMNS is formed, the shape of the 
spectrum is quite simple (see Fig. [6] (5b)). The effec- 
tive amplitude decreases monotonically until / ~ 5 kHz. 
Note that the power of the spectrum shape changes at 
/ ~ 3 kHz. This frequency is called the cut-off frequency, 
which is related to the compactness of two neutron stars 
(see Ref. [HI for details). The bump around / ~ 6 kHz 
is caused by the rotation of the merged object just before 
the collapse to a black hole (see Fig. [TJ top center). We 
find that this bump is enhanced in the case that the total 
mass is close to the critical mass, too — -Merit- 

Type II. A spectrum shape of type II is shown in Figs. [3] 
(2b) and[B](6b). In this case, we find several peaks caused 
by the rotation and oscillation modes of the HMNS. 
There is the maximum peak at / ~ 2.5-3 kHz caused 
by the rotation of the HMNS. The frequency of this peak 
is higher for the more compact HMNS, because the angu- 
lar velocity of the merged neutron stars is approximately 
written by oc m^ 2 / R^l 2 at the onset of the merger. 
Here R ns is the radius of two neutron stars. However, 
the peaks in the spectrum are too small to be detected 
by the advanced detectors such as advanced LIGO and 
LCGT for D = 100 Mpc, because the lifetime of the 
HMNS is short and the accumulated gravitational-waves 
cycles are small. There is a shallow decay of the spec- 
trum around 4 < / < 7 kHz. This is the common feature 
in the case that a HMNS is formed (see Figs.|S](2b)-(4b), 
and M (6b)). 
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Type III. Figures [5] (lb) and [5] (3b) show the spectrum 
shape of APR4-27 and H4-27. We find that the peak 
amplitude of the spectrum at / ~ 2.5 kHz is larger than 
that for type II and the primary oscillation mode appears 
clearly. Note that the amplitude of the peak is slightly 
smaller than the sensitivity curve of advanced LIGO (Op- 
timal NSNS version) and LCGT (Broadband version) for 
a hypothetical distance of 100 Mpc. If this peak could be 
observed with optimized detectors or for an event of short 
distance, we would get information about the physics 
of HMNSs and could constrain the supernuclear-density 
EOS H3- 

Note that the spectrum shape for type III depends on 
the EOS. For example, the spectrum shape of PS-27 is 
different from that for APR4-27 and H4-27 (see Fig. 
(lb), [5] (3b), and Fig. |5] (4b)). The bump around the 
maximum peak at / ~ 2.5 kHz in the spectrum of PS-27 
is broader than that of H4-27. 



IV. SUMMARY AND DISCUSSION 

We studied the dependence of the dynamical behav- 
ior of the binary neutron star merger on the EOS of 
the supernuclear-density matter in numerical relativity 
with piecewise-polytropic EOSs. We performed numer- 
ical simulations for 6 stiff EOSs and for several total 
masses in the range of 2.7M Q -3.0M Q . For all the cases, 
the equal-mass binary system was considered. 

We classified the merger process into three types: (i) 
a black hole is promptly formed (type I); (ii) a short- 
lived HMNS is formed (type II); (iii) a long-lived HMNS 
is formed (type III). For a given total mass of the bi- 
nary neutron star, the type of the merger process depends 
strongly on the EOS. In particular, the compactness of 
the neutron stars and the maximum mass of spherical 
neutron stars for a given EOS are key quantities for de- 
termining whether a HMNS is formed transiently or not. 
We found that the critical mass of the prompt formation 
of a black hole is in the range of 1.3 < M CT n/M max < 1.7. 
The latest observation of the high-mass neutron star with 
mass 1.97 ± O.O4M suggests that M crit is larger than 
2.7M . These facts indicate that a HMNS should be 
formed at least transiently for the total mass of the bi- 
nary neutron stars, tuq ~ 2.6M Q or less. We found that a 
long-lived HMNS is formed for tuq ~ 2.7 M@ and for the 
EOS with which the maximum mass of spherical neutron 
stars exceeds 2M Q , such as APR4 and H4. 

We studied the properties of a torus around a black 
hole formed after the merger. We found that the torus 
mass depends strongly on the type of the merger pro- 
cess. Specifically, the range of the torus mass is Mt orus < 
O.O1M for type I, 0.02M© < M torus < 0.05M© for type 
II, and 0.04 < M torus < 0.18M Q for type III. Thus we 
found that the torus mass is larger in the case that the 
HMNS is formed than in the case that a black hole is 
formed promptly. This is because materials in the outer 



envelope of the HMNS receive angular momentum from 
the central part via gravitational torques which work on 
the non-axisymmetric ellipsoidal HMNS. As a result, a 
part of the materials in the outer envelope does not fall 
into the black hole soon after the collapse of the HMNS 
and remains around the black hole to be a torus. In this 
sense, we conclude that the HMNS will play an important 
role for the merger scenario of short GRBs. Note that we 
performed numerical simulations only for the equal-mass 
system. For the unequal-mass case, the dynamical evo- 
lution of the merger could be modified. The heavier star 
may tidally disrupt the less massive companion if the 
mass ratio is high. As a result, the torus mass around 
the black hole may be enhanced (see Refs. (44T-|46j| ) . 

We also studied gravitational waves from the binary 
neutron star merger. The gravitational waveforms and 
their spectra depend strongly on the merger process. For 
type II and type III, we found that a large amount of 
gravitational waves are emitted by the HMNS which has 
a non-axisymmetric configuration. With decreasing its 
angular momentum, the HMNS approaches an axisym- 
metric ellipsoid and the amplitude of gravitational waves 
decreases. 

The amplitude of the spectrum for type I decreases 
monotonically with increasing / until 3-4 kHz because 
the HMNS is not formed. Note that there is a bump 
at / ~ 5-6 kHz caused by gravitational waves emit- 
ted by a merged object just before the collapse. By 
contrast, for type II and type III, the spectra have a 
complicated shape caused by the oscillation and rotation 
of the HMNS. At / ~ 2-3 kHz, there is the maximum 
peak cause by the fundamental rotation mode of the non- 
axisymmetric HMNS. We also found that several peaks 
due to the oscillation mode coupling with the rotation 
exist at both sides of the maximum peak. 

The detectability of gravitational waves from the 
HMNS is estimated. With advanced detectors such as 
advanced LIGO and LCGT, we may be able to detect 
gravitational waves caused by the rotation mode of the 
very long-lived HMNS (th > 10 ms) at around / ~ 2- 
3 kHz if the distance to the source is much smaller than 
100 Mpc. If these waves are observed, we will get in- 
formation about the physics of HMNSs and may be able 
to constrain properties of the supernuclear-density mat- 
ter go|. 



Acknowledgments 

We thank Y. Sekiguchi, Y. Suwa, and T. Naka- 
mura for useful discussions and comments. This work 
was supported by Grant-in-Aid for Scientific Research 
(21340051), by Grant-in-Aid for Scientific Research on 
Innovative Area (20105004), by the Grant-in-Aid of 
JSPS, by HPCI Strategic Program of Japanese MEXT, 
and by Grant-in-Aid for Young Scientists (B) 22740178. 



13 



[1] B. P. Abbott et al, Rep. Prog. Phys. 72, 076901 (2009). 
[2] T. Accadia et al., Class. Quantum Grav. 28, 025005 
(2011). 

[3] S. Hild, Class. Quantum Grav. 23, S643 (2006). 
[4] R. O'Shaughnessy, V. Kalogera, and K. Belczynski, As- 
trophys. J. 716, 615 (2010). 



[5] Adcanced LIGO, http://www.advancedligo.mit.edu/ [ 
[6] Advanced VIRGO, http://www.cascina.virgo.infn.it/advirgo/ 
[7] K. Kuroda et al, Class. Quantum Grav. 27, 084004 [ 
(2010). 

[8] R. Narayan, B. Paczynski, and T. Piran, Astrophys. J. | 

395, L83 (1992). 
[9] B. Zhang and P.Meazaros, Int. J. Mod. Phys. A 19, 2385 [ 

(2004). 

[10] T. Piran, Rev. Mod. Phys. 76, 1143 (2005). [ 

[11] S. Setiawan, M. Ruffert, and H.-Th. Janka, Mon. Not. [ 
R. Astron. Soc. 325, 753 (2004); W. H. Lee, E. Ramirez- 
Ruiz, and D. Page, Astrophys. J. 632, 421 (2005); M. [ 
Shibata, Y. Sekiguchi, and R. Takahashi, Prog. Theor. 
Phys. 118, 257 (2007). [ 

[12] M. Shibata, K. Taniguchi, and K. Uryu, Phys. Rev. D 

68, 084020 (2003). [ 

[13] M. Shibata, K. Taniguchi, and K. Uryu, Phys. Rev. D 

71, 084021 (2005). [ 

[14] M. Shibata and K. Taniguchi, Phys. Rev. D 73, 064027 
(2006). 

[15] L. Biotti, B. Giacomazzo, and L. Rezzolla, Phys. Rev. D [39 
78, 084033 (2008). [40 

[16] K. Kiuchi, Y. Sekiguchi, M. Shibata, and K. Taniguchi, [41 
Phys. Rev. D 80, 064037 (2009). [42 

[17] R. Oechslin, H.-Th. Janka, and A. Marek, Astorn. As- 
trophys. 467, 395 (2007). 

[18] J. S. Read, B .D. Lackey, B. J. Owen, and J. L. Friedman, [43 
Phys. Rev. D 79, 124032 (2009). 

[19] F. Ozel and D. Psaltis, Phys. Rev. D 80, 103003 (2009). [44 

[20] T. Yamamoto, M. Shibata, and K. Taniguchi, Phys. Rev. 

D 78, 064054 (2008). [45 

[21] M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 

(1995). [46 

[22] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 

024007 (1998). [47 

[23] M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlo- 
chower, Phys. Rev. Lett. 96, 111101 (2006); J. G. Barker, 
J. Centrella, D.-I. Choi, M. Koppitz, and J. van Meter, [ 
Phys. Rev. Lett. 96, 111102 (2006); L. Baiotti, M. Shi- 
bata, and T. Yamamoto Phys. Rev. D 82, 064015 (2010). 

[24] B. Briigmann, J. A. Gonzalez, M. Hannam, S. Husa, U. [49 
Sperhake, and W. Tichy, Phys. Rev. D 77, 024027 (2008). 

[25] A. Kurganov and E. Tadmor, J. Comput. Phys. 160, 241 [50 
(2000). 



A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, 
Phys. Rev. C 58, 1804 (1998). 

F. Douchin and P. Haensel, Astron. Astrophys. 380, 151 
(2001). 

B. D. Lackey, M. Nayyar, and B. J. Owen, Phys. Rev. D 
73, 024021 (2006). 

N. K. Glendenning and S. A. Moszkowski, Phys. Rev. 
Lett. 67, 2414 (1991). 

M. Alford, M. Brady, M. Paris and S. Reddy, Astrophys. 
J. 629, 969 (2005). 

V. R. Pandharipande and R. A. Smith, Nucl. Phys. 
A237, 506 (1975). 

P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. 
Roberts, and J. W. T. Hessels, Nature 467, 1081 (2010). 
J. M. Lattimer and M. Prakash. larXiv: 1012.32081 1 (2010). 
F. Ozel, D. Psaltis, A. Ransom, P. Demorest, and M. 
Alford, Astrophys. J. 724, L199 (2010). 
A. Bauswein, H.-Th. Janka, and R. Oechslin, Phys. Rev. 
D 82, 084043 (2010). 

Y. Sekiguchi, K. Kiuchi, K. Kyutoku, and M. Shibata, 
submitted to Phys. Rev. Lett. (2011). 
K. Taniguchi and M. Shibata, Astrophys. J. Suppl. 188, 
187 (2010). 

E. Gourgoulhon, P. Grandclement, K. Taniguchi, J. -A. 
Marck, and S. Bonazzola, Phys. Rev. D 63, 064029 
(2001). 

L. Bildsten, and C. Cutler, Astrophys. J. 400, 175 (1992). 

C. S. Kochanek, As trophys. J. 398, 234 (1992). 
LORENE webpage: http://www.lorene.obspm.fr/ 

C. M. Zhang, J. Wang, Y. H. Zhao, H. X. Yin, L. M. 
Song, D. P. Menezes, D. T. Wickramasinghe, L. Ferrario, 
and P. Chardonnet, Astron. Astrophys. 527, A83 (2011). 
T. W. Baumargte, S. L. Shapiro, and M. Shibata, Astro- 
phys. J. 528, L28 (2000). 

K. Kiuchi, Y. Sekiguchi, M. Shibata, and K. Taniguchi, 
Phys. Rev. Lett. 104, 141101 (2010). 
M. Shibata, K. Taniguchi, and K. Uryu, Phys. Rev. D 
68, 084020 (2003). 

L. Rezzolla, L Baiotti, B Giacomazzo, D. Link, and J. A. 
Font, Class. Quantum Grav. 27, 114105 (2010). 
M. Boyle, D. A. Brown, L. E. Kidder, A. H. Mroue, H. P. 
Pfeiffer, M. A. Scheel, G. B. Cook, and S. A. Teukolsky, 
Phys. Rev. D 76, 124038 (2007). 

For planned sensitivity curves for ad- 
vanced LIGO, see https://dcc.ligo.org/cgi- 
bin/DocDB/ShowDocument?docid=2974 
M. Punturo et al. Class. Quantum Grav. 27, 194002 
(2010). 

M. Shibata, Phys. Rev. Lett. 94, 201101 (2005). 



